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Abstract 

In the single degenerate (SD) scenario of type la supernovae (SNe la), the collision of the 
ejecta with its companion results in stripping hydrogen rich matter from the companion star. 
This hydrogen rich matter might leave its trace in the light curves and/or spectra. In this paper, 
we perform radiation hydrodynamical simulations of this collision for three binary systems. As 
a result, we find that the emission from the shock-heated region is not as strong as in the 
previous study. This weak emission, however, may be a result of our underestimate of the 
coupling between the gas and radiation in the shock interaction. Therefore, though our results 
suggest that the observed early light curves of SNe la can not rule out binary systems with a 
short separation as the progenitor system, more elaborate numerical studies will be needed 
to reach a fair conclusion. Alternatively, our results indicate that the feature observed in the 
early phase of a recent type la SN 2014J might result from interaction of the ejecta with a red 
giant companion star. We also discuss the dependence of spectral features of Ha and Si II 
absorption lines on viewing angles and investigate whether they can constrain the event rate 
of the SD progenitor. 


1 



Key words: hydrodynamics - line: profiles - radiative transfer - supernovae: general 


1 Introduction 

Type la supernova (SN la) is thought to originate from an accreting carbon-oxygen white dwarf (WD) 
in a binary system (Whelan & Iben 1973; Iben & Tutukov 1984; Nomoto et al. 1984; Webbink 1984). 
However, it is uncertain what is the companion star supplying its mass. There are two scenarios 
to explain SNe la. One is the single degenerate (SD) scenario, in which the companion star is a 
main sequence or a red giant star (Iben & Tutukov 1984; Nomoto et al. 1984). The envelope of 
the companion star overflows the Roche lobe, and transfers to the surface of the WD. The other is 
the double degenerate (DD) scenario, in which two WDs merge and end up with explosion (Iben & 
Tutukov 1984). To distinguish the two scenarios, we need to find observable differences between 
them. 

In the SD scenario, SN ejecta collide with its companion. Kasen (2010) simulated this col¬ 
lision and found that high energy photons are emitted from the shock-heated region into a certain 
limited solid angle and that the emission becomes prominent before the peak of the light curve due 
to 56 Co decay. This emission is strong especially in a WD + red giant system. Hayden et al. (2010) 
investigated early phases of observed light curves gathered by the Sloan Digital Sky Survey (SDSS) to 
see whether they have this prompt emission. By comparing its early light curve with the observations, 
they found that there is no obvious emission feature in the light curves, which constrained the event 
rate of each model. As a result, a progenitor system with a main sequence companion more massive 
than 6 M 0 or a red giant is ruled out as the primary source of SNe la. 

Kasen (2010) assumed the LTE condition and further that whole matter is composed of radi¬ 
ation dominated gases with 7 = 4/3. We suspect that they overestimated the radiation temperature 
of the emission regions. In fact, based on their detailed calculations of spectra from the same hydro- 
dynamical models presented here, Maeda et al. (2014) argued that hydrogen rich matter filling in the 
hole excavated by the companion star prevents the photosphere from quickly receding to the Ni-rich 
region and makes the SN redder than expected by the previous study (Kasen et al. 2004). In addition, 
Kasen (2010) focused only on the impacts on light curves. We investigate also the influence on spec¬ 
tra. If the material stripped from the companion does not spread in a wide solid angle, hydrogen lines 
from the companion might appear in the spectrum, depending on the line of sight. 
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Leonard (2007) investigated nebula spectra of some SNe la and found the upper limit 
of 0.01 M 0 for solar abundance material. Several authors (Wheeler et al. 1975; Marietta et al. 
2000; Pakmor et al. 2008) calculated the mass of the material stripped from the companion star by 
the collision. Marietta et al. (2000) performed hydrodynamical simulations to estimate the stripped 
mass. According to their results, the stripped masses are 0.15M e for the 1M 0 main sequence com¬ 
panion, O.53M 0 for the 1M 0 red giant companion. On the other hand, there is no spectral feature 
of the stripped material in observed spectra. Recently, Pakmor et al. (2008) adopted a more realistic 
profile for a main sequence companion from a binary evolution theory and found that the stripped 
mass decreases to 0.01 — O.O2M 0 . (Kasen et al. 2004) investigated signatures of hydrogen lines in 
the phase later than the maximum light when the shock heating due to the collision has been already 
consumed by adiabatic expansion. 

In this study, we relax the assumption of the instantaneous coupling between gas and radiation 
to examine whether or not the decoupling of gas and radiation changes the influence on light curves 
found by Kasen (2010). In addition, we investigate whether hydrogen in the stripped material leaves 
its trace in spectra especially before maximum light. In later phases, Maeda et al. (2014) already 
studied spectra for the same model presented here. At first, we simulated the collision of a SN la with 
its companion with a radiation hydrodynamical code. Secondly, we calculated spectra by ray-tracing 
snapshots of the distributions of density, temperatures, and velocity obtained from the simulation. 
Section 2 describes our models. Section 3 and 4 present the numerical methods for the simulation and 
ray-tracing. In section 5, we describe our results. Section 6 concludes this study. 

2 Models 

As progenitor systems, we consider three models named MS, RGa, and RGb (Table. 1). Model MS is 
a close binary system with a main sequence companion. We use a polytropic star with the index n = 3 
to model the structure of the companion star. We assume the companion star composed of a solar 
abundance material. Its mass M 2 is equal to 1M 0 . The binary separation A is equal to 3 x 10 11 cm. 
This model is the same as the 1 M & main sequence model in Kasen (2010). 

Model RGa is a close binary system with a red giant companion. In this model, the companion 
star consists of the helium core and the hydrogen-rich envelope with the solar abundance. Its core 
mass and total mass M 2 are equal to O.4M 0 and 1M 0 , respectively. The binary separation A is 
equal to 2 x 10 13 cm. We assume a fully convective envelope with n = 3/2. In the hydrodynamical 
simulation, the core is treated as a point source of gravity. This model is the same as the 1M 0 red 
giant model in Kasen (2010). Model RGb is the same as model RGa but with a longer separation, i.e., 
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Table 1 . Models in our calculation. 


Model 

M 2 (M q ) 

R .2 (10 13 cm) 

A (10 13 cm) 

MS 

1 

0.01 

0.03 

RGa 

1 

0.7 

2 

RGb 

1 

1 

3 


A = 3 X 10 13 cm. 

The radii R 2 of the companion stars in the above models are obtained by an empirical law for 
close binary systems (Paczynski 1971) assuming that the companion stars fill the Roche lobe. 

~r = 0-38 + 0-2 log ( 1 ) 

A lvl wd 

Here, Mwd is the mass of the progenitor WD, which is equal to 1.38M 0 . Thus R 2 = 0.35A 


3 Simulation 


We simulate the collision between a SN la and its companion. This simulation is performed with 
a radiation hydrodynamical code in the axisymmetric spherical coordinate system (r, 9). We ignore 
the orbital motion of the progenitor binary. The orbital velocities are calculated as 190 km/sec for 
model MS, 23 km/sec for model RGa, and 19 km/sec for model RGb. The material expanding much 
faster than these velocities is not affected by the orbital motion. Since the material slower than these 
velocities is located inside the photosphere over the relevant period, we do not observe the influence 
of the orbital motion. 

The basic equations governing hydrodynamics in an Eulerian coordinate system are written as 
follows. 


dp 

~di 

d_ 

dt 

d_ 

dt 


+ V • (pv) = 0, 


(pv) + V • (pvv) = 



+ V- 


-Vp + fe x t, 



v ' fext + <5, 


( 2 ) 

( 3 ) 

( 4 ) 


where t is the time, p the density, v the fluid velocity, p the gas pressure, e the specific internal energy 
of gas, h = e+p the specific enthalpy, f ext the external force, and S is the energy generation rate from 
heat sources. For ideal gases, p and e are expressed by the equation of state with the adiabatic index 


7, 

pkT 

P= -, 

pm p 


(5) 
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(6) 


7 pkT 

e= - 7 -’ 

7-1 jj/rn p 

where k is the Boltzmann constant, m p the proton mass, // the mean molecular weight, and T is the 
temperature. In addition, the density of each isotope can be calculated by equations 

^(A»+V-(A>v) = A, (7) 


where X\ and I), are the mass abundance and the production/decay rate of the z-th isotope. 
Simultaneously, we solve the bolometric radiation energy equation, assuming the flux limited dif¬ 
fusion for the radiation flux F. 


dE 

~dt 

F = 


V • [F + AE ) = -v • -F + ca(aT 4 - E), 
-cWE 


3y+| VE\/E* 


( 8 ) 

(9) 


where c is the speed of light in vacuum, a the radiation constant, x — Oi + a,a the absorption co¬ 
efficient, <7 the scattering coefficient, and E is the radiation energy density. The flux limited diffu¬ 
sion approximation has been widely used in multi-dimensional radiation-hydrodynamics simulations 
(Ohsuga et al. 2009, for example). Considering the interaction between gas and radiation, f ext and S 
are written as, 


fext = pg + -F, (10) 

c 

S = ca(E — aT 4 ) + e, (11) 


where g is the gravitational field, and e is the heating rate by nuclear reactions. 

To solve the above system of equations, we separately integrate each term with respect to time 
according to the operator splitting method with first order accuracy in time. We integrate the advection 
terms with the Harten-Lax-van Leer-Contact (HLLC) scheme (Toro et al. 1994). In addition, we 
use the Monotonic Upstream-Centered Scheme for Conservation Laws (MUSCL) scheme (van Leer 
1976) to obtain solutions with the second order accuracy in space. We solve the radiative diffusion 
and the thermalization processes with an implicit method. 


3.1 Initial condition 

A companion star is located at the origin of the coordinate system. We put the center of the expanding 
ejecta at a point r = A, 6 = 0. The initial distributions of the density and abundance of each element in 
the expanding ejecta are obtained from W7 model (Nomoto et al. 1984). They are shown in Figure 1. 
Considered isotopes are 12 C, 16 0, 20 Ne, 24 Mg, 28 Si, 32 S, 36 Ar, 40 Ca, 52 Fe and 56 Ni. Calculations are 
initiated at the moment when the ejecta touch the surface of the companion. Figure 2 shows the initial 
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density distribution of each model. 


3.2 Input physics 

Here, we describe how to calculate gravity, transport coefficients, and nuclear reactions in equa¬ 
tions (2)-(ll). Since the gravity mainly affects the amount of the stripped material, we assume that 
the companion star is the only source of gravity, and further that the gravitational field does not change 
with time. About the transport process, we take into account the free-free absorption and the Thomson 
scattering. We use the Planck-mean of the free-free absorption coefficient as the bolometric absorp¬ 
tion coefficient. It is underestimated because we ignore absorption by bound-bound and bound-free 
transitions. 

The ejecta are heated by gamma-ray from the radioactive decays of 56 Ni to 56 Co, and 56 Co to 
56 Fe. The life time of 56 Ni we use in our simulations is r( 56 Ni) = 8 days and the energy emitted from 
its decay 77( 56 Ni) = 2.13 MeV (Amett 1996). For the 56 Co decay, we use r( 56 Co) = 111 days and 
E( 56 Co) = 4.56 MeV. The decay rates of isotopes are written as 


D( 56 Ni) = -m( 5b Ni) 


56 


n 


( 56 Ni) 


56ru\_ n ( Ni) /56 n~\ n ( ^°) 


D( 5b Co) = m( 5 b Ni) 
D( 56 Fe) = m( 56 Co) 


r( 56 Ni) ’ 

56 Mn n f56 

r( 56 Ni) _m ( 56 C °)r( 56 Co) ?: 
n( 56 Co) 


( 12 ) 

(13) 

(14) 


r( 56 Co) ’ 

where m(A) and n(A) are the mass and number density of the isotope A, respectively. Then, the 
gamma-ray energy generation rate S 1 is calculated by 


£,=B("Ni)S+®rc»)" ( “ 0o) 


(15) 


r( 56 Ni) v r( 56 Co) ’ 

We do not calculate the transport of emitted gamma-ray photons. Instead, we assume that a fraction 
of the generated gamma-ray is immediately absorbed by the local material to the extent determined 
by the optical depth of the ejecta. The heating rate e is then written as 


e = £ 7 [1 — exp(—r 7 ) . 


(16) 


Here r 7 is the gamma-ray optical depth calculated by 

r 7 = J Kypdr, (17) 

where k 7 is the gamma-ray absorption opacity. Here the optical depth is integrated along the line 
passing the center of the ejecta. Fransson & Chevalier (1989) studied the gamma-ray transport in 
SN la, assuming the gamma-ray is not scattered, only absorbed by the ejecta. They obtain k 7 = 
0.03 cm 2 g -1 to reproduce the observed light curves of SNe la. Although our assumption might be 
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a crude approximation, the total amount of thermal energy converted from the gamma-ray energy is 
almost correct in the sense that it can reproduce the typical light curve of SN la and the result of 
Maeda et al. (2014). 

3.3 Remapping 

Initially, the physical scale of the binary system is 10 11 cm in model MS, and 10 13 cm in models RGa 
and RGb. On the other hand, at the maximum light (~ 20 days after the explosion), the radius of the 
ejecta becomes 3 x 10 15 cm. In order to trace the profiles of the ejecta and the stripped material, we 
expand the calculation region according to the expansion of the ejecta. When the size of the ejecta 
reaches 90% of the size of the calculation region, we increase the r-axis interval of computational 
grids by a factor of two. The values on the outermost cells of the original region are substituted to 
fluid variables in newly created cells outside the original calculation region. 


4 Ray tracing 

As results of the above simulation, we obtain snapshots of the spacial distribution of each variable. 
Integrating the radiation transport equation along a specified line of sight on the snapshot, we obtain 
light curves and spectra. Here, we define the Cartesian coordinates (x, y , z) of which the 2 -axis 
is along the line of sight. Considering both thermal emission and isotropic scattering, the surface 
brightness along the 2 -axis is written as 

I„(x,y) = J{a u + cr u )(S u - I v )dz, (18) 

where is the intensity along the 2 -axis. S v is the source function defined by 

a u B u (T) + cr u J u 


S v = 


(19) 


CYjj <Jv 

where B U (T) is the Planck function, and ./„ is the mean intensity. Then, the isotropic luminosity is 
calculated by 


LT = 4tt 


l v (x,y) dxdy. 


( 20 ) 


From the radiation hydrodynamic simulation, we have distributions of the bolometric mean 
intensity J = (c/Att)E. Then, the bolometric source function S is calculated by 

acaT^/An + aJ 


S = 


( 21 ) 


a + cr 

Here, a and a the bolometric absorption and scattering coefficients, which are equal to those used 
in the above simulation. Computing the bolometric luminosity from each snapshot, we obtain the 
bolometric light curve. 
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Table 2. Parameters of spectral lines. — / is the oscillator strength. g. t and E t are the statistical weight and the level energy of the 


lower level i. gj and Ej are those of the upper level j. 


Lines 

'og (gf) 

Ei (eV) 

Eg (eV) 

9i 

9j 

HI A6562.81A 

0.710 

10.198 

12.087 

8 

18 

C II A6578.05A 

-0.021 

14.449 

16.333 

2 

4 

C II A6582.88A 

-0.323 

14.449 

16.332 

2 

2 

Sill A6347.10A 

0.149 

8.121 

10.073 

2 

4 

Si II A6371.36A 

-0.082 

8.121 

10.066 

2 

2 


To obtain a spectrum, we need to calculate the monochromatic source function S v . Though 
we need the monochromatic mean intensity J u to calculate S v , we know only the value of bolometric 
mean intensity J from the simulation. To evaluate J u , we assume that J u takes a form of the Planck 
function, 

J u = B u (T iad ), J = ca.T r 4 ad /4vr. (22) 

To determine the absorption and scattering coefficients, we take into account Thomson scattering, 
free-fee absorption, and bound-bound absorption by relevant atomic lines. Ionization states of each 
isotope are calculated by the Saha equation. To calculate the optical depth of the line absorption in 
moving media, we apply the Sobolev approximation (Rybicki & Hummer 1978). In this study, we 
only concern with hydrogen features in optical spectra. Then, we calculate only Ha, C II, and Si II 
lines. These lines are located in a narrow range of the rest wavelengths between 6300 A and 6600 A. 
Table 2 shows spectral lines considered in this study. 

When neutral hydrogen is located outside the photosphere where the temperature is lower than 
the photospheric temperature, Ha line might appear as absorption in the spectrum. 

5 Results 

Results from the simulation and ray-tracing explained in the previous sections are presented in the 
following subsections. 

5.1 Stripped mass 

From the above simulation, we obtain spatial structures of the ejecta and the stripped material. The 
stripped mass from the companion star is calculated by adding up the mass of hydrogen rich gas with 
the fluid velocity exceeding the escape velocity. The masses are O.45M 0 in model MS, O.51M 0 in 
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model RGa, and O.49M 0 in model RGb. In model MS, the amount is overestimated, compared to 
the previous studies (about 0.01 — O.O2M 0 in Pakmor et al. (2008)). This is because we do not have 
enough resolution in the low velocity region around the companion star due to the remapping process. 
On the other hand, in models RGa and RGb, the amount is slightly smaller than the previous studies 
(Wheeler et al. 1975; Marietta et al. 2000). Nevertheless, the trend that the entire envelope is stripped 
by the impact is consistent with these studies. 


5.2 Structures obtained from the simulation 


Figures 3-5 show distributions of density p and total number abundance X/(2Q/A) at days 1 and 20 
after the explosion in each model. As shown in the figures, the hydrogen-rich material is stripped 
along the line 6 = n. At day 20, the ejecta and the stripped material are already expanding homolo¬ 
gously. A cavity is created around the axis 9 = it with the spread angle of ~ 7 t/6 rad. If we view this 
SN with a viewing angle in this range, the stripped hydrogen is not hidden by the surrounding ejecta. 

Figure 6 shows one-dimensional structures of the density p, the gas temperature T, and the 
radiation temperature T rad along the line 0 = tt at day 1. We calculate the position of the photosphere 
where thermalization length is equal to unity, which is indicated by the dotted vertical line. The 
photospheric temperature is equal to 1.0 x 10 4 K for model MS, 5.2 x 10 4 K for model RGa, and 

6.2 x 10 4 K for model RGb. Kasen (2010) analytically estimated the effective temperature of the 
prompt emission as 


T eS = 2.5 x 10 4 a{i 4 K- 35 / 3 %£ /72 K, 


(23) 


where a, 3 = A/10 13 cm, n e is the electron scattering opacity, and t day is the time in days. At t = 1 day, 
T eS is estimated as 5.0 x 10 4 K for model MS, 1.4 x 10 5 K for model RGa, and 1.6 x 10 5 K for 
model RGb, In comparison, the temperatures become lower in our simulation. This is because the 
collision heats the gas first in our simulations, which separately treat the temporal evolutions of the 
radiation and gas energy densities. The subsequent emission from the gas is not enough to raise the 
photospheric temperature as high as this formula due to the low densities. 


5.3 Light curves 

Bolometric light curves with viewing angles 6 = 0, ir/2, 2 tt/3, 57t/6 and tt are shown in Figure 7. 
The collision radiates prompt emission especially in the range of 6* > 5zr/6. Kasen (2010) analytically 
estimated bolometric luminosities due to the collision at day 1 to be 3 x 10 41 erg s _1 for the 1M 0 
main sequence model, 2 x 10 42 erg s _1 for the 6M 0 main sequence model, and 2 x 10 43 erg s _1 for the 
1 M q red giant model. In comparison with our results, the luminosity is equal to 6.9 x 10 40 erg s _1 for 
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model MS, 2.3 x 10 42 erg s -1 for model RGa, and 3.4 x 10 42 erg s^ 1 for model RGb. The luminosities 
in models MS, RGa are dimmer than the previous study. This is because the temperature of emission 
region becomes lower than their results. In the analysis of Hayden et al. (2010), the luminosity in 
the 6M 0 main sequence model is the threshold above which a progenitor system is ruled out from 
the primary source of SNe la. The luminosity in model RGa is on this threshold. A model with a 
slightly longer separation as model RGb is to be ruled out because the emission can yet be detected. 
Here one should note that our models may underestimate the strength of the coupling between the 
gas and radiation because the free-free emission and the energy exchange by Compton scattering, 
which our models do not include, accelerate the equilibration as pointed out by Weaver (1976); Katz 
et al. (2010). As a result, we may significantly underestimate the luminosity due to the shock heating. 
More elaborate numerical studies including the Comptonization are needed to assess this effect. In 
Figure 8, we compare our results of model light curves in the R and V bands viewing from the side of 
the companion star (6 = ^ r) with two well observed type la SNe 201 lfe and 2014J. Though the early 
light curve of SN 2014J is of unfiltered observations or of approximate R magnitudes estimated from 
narrowband detections using H Q filters(Goobar et al. 2014a; Goobar et al. 2014b), the shape indicates 
some influence from the collision between the SN ejecta and the companion star. Therefore SN 2014J 
might have a red giant as the ex-companion star and prefers the SD scenario. Note that we do not 
intend to fit the observed light curve of SN 2014J with our current models. A detailed comparison 
will be discussed elsewhere. 

5.4 Spectra 

We calculate spectra with several viewing angles 9 every 5 days up to 40 days since explosion. 
Figures 9-11 show spectrum of each model. Here the monochromatic luminosity L\ is normalized 
by the continuum luminosity L c x . In the early spectra, strong absorption features come from C II line 
(~ 6300 A) and Si II line (~ 6100 A). On the other hand, Ha feature is quite weak or hidden by the C 
II absorption, because W7 model retains a large amount of carbon in the outer ejecta. In reality, C II 
absorption features are weak and rarely detected in observed SN la spectra (Folatelli et al. 2012). We 
overestimate their absorption feature and therefore there is a possibility to detect Ha feature in some 
SNe la depending on the carbon content and the viewing angle. 

Especially, model MS has a strong Ha feature on its spectra before day 10 and after day 30. 
As was mentioned before, the overestimate of the stripped mass in this model enhances the feature as 
compared with reality. In models RGa and RGb, Ha absorption becomes prominent only after day 
30. On the other hand, C II feature already becomes weak in this epoch and does not disturb the Ha 
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feature. The S/N ratio needed for the 3a detection of these features is written as 

S/N = 3x (24) 

\ L * ~L X \ 

At day 35, the S/N ratio to detect the most strong feature in the spectrum is calculated as 59 for model 
RGa, and 125 for model RGb. Figure 12 shows the time evolution of the equivalent width of Ha line. 
The peak value of the equivalent width with the viewing angle 6 = 7r is 6 A for model MS, 5 A for 
model RGa, and 4 A for model RGb. With 6 = 57t/6, the peak value decreases by a factor of about two 
compared with the above value in each model. Thus, in models RGa and RGb, the Ha feature can be 
detected when the viewing angle is in the range of 6 > II-k/YI, or within the solid angle of H ~ 0.2. 
The event ratio in each model is estimated as Vt/An ~ 0.02. However, this Ha feature has never been 
observed in SN la spectra. Thus, these systems are ruled out from the progenitor candidates. Here, 
it should be noted that these calculations assume that the ionization states are in thermal equilibrium 
and that we need to know the continuum spectra with sufficient accuracies in advance. Maeda et al. 
(2014) already discussed hydrogen lines in this late phase in more detail including the Paschen lines. 

As another spectral feature to be investigated, Figure 13 shows time evolution of silicon ve¬ 
locity with different viewing angles. The differences of velocities between 6 = 0 and n are about 
1000 km s -1 for all models. A spectral resolution of about 20 A could distinguish the differences. 
This variation occurs because the edge of the ejecta is decelerated and masked by the companion star. 
We usually observe a SN from a single direction and it is difficult to distinguish this effect from the 
variety of expansion speeds of individual SNe la. One possibility to extract this effect is a SN la for 
which a several light echoes are discovered. They will have different silicon velocities depending on 
the viewing angles, and can be compared with our results. Krause et al. (2008) obtained a light echo 
spectrum of SNR Tycho with the spectral resolution of 24 A, observed by the FOCAS camera on the 
Subaru telescope. The resolution we require is higher. Furthermore, the light echo spectrum is con¬ 
structed by a mixture of SN spectra for about 10 days. Thus, it is not feasible to detect the variation 
at present. 

These spectral features confirmed the results obtained by more sophisticated treatment of ra¬ 
diative transfer in the same models (Maeda et al. 2014) except that the temporal evolution of sili¬ 
con line. Maeda et al. (2014) obtained the silicon velocities independent of viewing angles at day 
~ 10 while we always see lower velocities when viewing from the companion side. This difference 
might be partly due to different treatments of the shock heating. Maeda et al. (2014) redetermined 
the temperature by imposing the radiative equilibrium condition ignoring the shock heating. Higher 
temperature due to the shock heating reduces the number of Si II ions and lowers the velocities at 
the absorption minima. In later phases, the adiabatic cooling diminishes contributions of the shock 
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heating and the results of the two different simulations converge. 


6 Summary and Conclusions 

In the SD scenario of SNe la, the collision between the ejecta and its companion occurs. Due to this 
collision, their light curve and spectra are changed from the case of SN la without the companion. 
These characteristics are expected to constrain progenitor systems. In this study, we perform radiation 
hydrodynamical simulations of this collision for three binary systems in the SD scenario. Using the 
obtained snapshots, we solve the radiation transport equation, and obtain their light curves and spectra 
with several viewing angles. Then, we determine whether they can constrain the event rate of the SD 
progenitor. The results are as follows. 

In comparison to Kasen (2010), the prompt emission due to the collision becomes weaker, 
since the photospheric temperature becomes lower by the weaker coupling of gas and radiation. In 
models MS and RGa, the luminosity does not exceed that of the 6M 0 main sequence model in Kasen 
(2010), above which the emission can be detected in observed SNe la. On the contrary, in model RGb, 
the luminosity still exceeds the threshold. Thus, we conclude that binary systems with the separation 
A < 2 x 10 13 cm is allowed as a progenitor of SNe la. As exemplified by SN 2014J, photometric 
observations during the first few days for a large sample of SNe are crucial for distinguishing the two 
scenarios. 

In all models, the equivalent widths of Ha absorption become a few A 30 days after the 
explosion. However, in model MS, the stripped mass is overestimated and the absorption will become 
weaker in reality. On the other hand, in models RGa and RGb, the Ha feature can be detected if the 
viewing angle 6 > ll7r/12. Then, 2% of SNe la in these models should have the feature on their 
spectra, while there is no sign of Ha absorption in the observed spectra. If a sufficient number of 
spectra were obtained with good S/N ratios and were not showed this feature, a model with a red 
giant companion filling the Roche lobe could be ruled out from SNe la progenitors. 

The central wavelength of Si II absorption varies depending on viewing angles. If more than 
one light echoes are discovered in a single SN, we obtain its spectra with different viewing angles. 
Since the difference of silicon velocities depending on the viewing angles is about 1000 km s _1 , the 
required resolution is about 20 A. The current observation has not reached this resolution for light 
echo spectra. We need higher resolution which will be reached by future observations. 

As shown above, our models indicate that the prompt emission can not constrain binary sys¬ 
tems with a short separation. One should note here that our models may significantly underestimate 
the luminosity due to the shock heating because we do not take into account the coupling of radia- 
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tion and gas through Compton scattering. Thus we need to address the effects of Compton scattering 
before reaching a conclusion on this issue. At the same time, we also found that the spectral fea¬ 
ture might constrain the occurrence of thermalization through Compton scattering. Furthermore, Ha 
feature brings us additional information of the progenitor of an individual SN la. Because we can 
not predict where and when a SN will occur, the prompt emission is difficult to be found. It will be 
detected only when a monitoring telescope is occasionally viewing around its site. In contrast, some 
hydrogen features become prominent after the maximum light, which relaxes the required observing 
cadence. If we continuously observe the SN and take spectra with a right viewing angle, we will never 
miss it. This is a great advantage. 
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Fig. 1. Density and abundance structure of the ejecta at 0.1 days after the explosion. 



Fig. 2. Initial condition of each model. — The left panel shows model MS. The center panel shows model RGa. The right panel shows model RGb. Note that 
scales are different in these models. 



Fig. 3. Density and abundance structures in model MS at days 1 and 20 after the explosion. — In each panel, the number abundance is shown on the left 
side, and the density is on the right side. The red region on the left side is rich in hydrogen. 
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Fig. 4. Density and abundance structures in model RGa at days 1 and 20 after the explosion. — In each panel, the number abundance is shown on the left 
side, and the density is no the right side. The red region on the left side is rich in hydrogen. 



Fig. 5. Density and abundance structures in model RGb at days 1 and 20 after the explosion. — In each panel, the number abundance is shown on the left 
side, and the density is on the right side. The red region on the left side is rich in hydrogen. 
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Fig. 6. Density and temperature structures along 6 = tt at day 1 after the explosion in each model. — The top panel shows model MS. The center panel 
shows model RGa. The bottom panel shows model RGb. The solid line is density, long-dashed line is gas temperature, and short-dashed line is radiation 
temperature, respectively. The dotted line shows position of the photosphere. 
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Fig. 7. Bolometric light curves with different viewing angles. — The top panel shows model MS. The center panel shows model RGa. The bottom panel 
shows model RGb. The dash-dotted, dotted, short-dashed, long-dashed, and solid lines show light curves when we observe the SN with the viewing angle 
9 = 0 , 7t/2, 2-7t/3, 5ir/6, and rr, respectively. 
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Fig. 8. Model light curves in the R band (top panel) and V band (bottom panel) viewing from the side of the companion stars. Light curves of well observed 
type la supernovae 2011 fe (Richmond & Smith 2012) and 2014J (Zheng et al. 2014) are also plotted for comparison. Note that the points of SN 2014J are of 
unfiltered observations except for filled squares in the top panel. The filled squares denote the R magnitudes estimated from detections using the narrowband 
He filters. The distance to SN 2014J is assumed to be 3.5 Mpc. 
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Fig. 9. Spectra normalized by the continuum for model MS at days 5, 10, 15, 20, 25, 30, 35 and 40 after the explosion. — The black, red, green, blue, and 


magenta lines show spectra when we observe the SN with viewing angles 8 = 0, 2rr/3, 5rr/6, lln-/l2, and -rr, respectively. The solid lines are calculated 


including all spectral lines mentioned in the text, while the dashed lines without Ha lines. 
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Fig. 10. Spectra in model RGa at day 5, 10, 15, 20, 25, 30, 35 and 40 after the explosion. — The black, red, green, blue, and magenta lines show light curves 
when we observe the SN with the viewing angle 9 = 0, 2tt/3, 5v/ 6, lln/12, and v, respectively. The solid lines are calculated including all spectral lines. 
Besides, the dashed lines are calculated without Ha lines. They are normalized by the continuum spectra. 
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Fig. 11 . Spectra in model RGb at day 5,10,15, 20, 25, 30, 35 and 40 after the explosion. — The black, red, green, blue, and magenta lines show light curves 


when we observe the SN with the viewing angle 9 — 0, 2tv/3, 5v/ 6, lln/12, and rr, respectively. The solid lines are calculated including all spectral lines. 


Besides, the dashed lines are calculated without Ha lines. They are normalized by the continuum spectra. 
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Fig. 12. Equivalent width of Ha line in each model. — The top panel shows model MS. The center panel shows model RGa. The bottom panel shows model 
RGb. The dotted, short-dashed, long-dashed, and solid lines show light curves when we observe the SN with the viewing angle 9 = 27t/3, 5-k/ 6, ll-ir/ 12 , 
and 7 r, respectively. 
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Fig. 13. Angular dependence on Si li line in each model. — The top panel shows model MS. The center panel shows model RGa. The bottom panel shows 
model RGb. The dash-dotted, dotted, short-dashed, long-dashed, and solid lines show light curves when we observe the SN with the viewing angle 0 = 0, 
7t/2, 27t/ 3, 5rr/6, and tt, respectively. 
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